Superconductivity and Topological Numbers in 
2^ ! the Hubbard Chain with Bond-Charge Interaction 

a^ 

A.A. Aligia, K. Hallberg, CD. Batista and G. Ortiz* 

O . 

^^ I Centra Atomico Bariloche and Instituto Balseiro, Comision Nacional de Energia 

Atomica, 84OO S.C. de Bariloche, Argentina 
* Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 8754-5 



(N 



> 

O 
00 



0^ 
0\ 



We determine the quantum phase diagram of the Hubbard chain with 
electron-hole symmetric correlated hopping at 1/2- and 1/4-filling using geo- 
metric concepts and continuum limit field theory. The long distance behavior 
of various correlation functions show a very rich phase diagram with several 

2 ! insulating, metallic, and superconducting phases, which might be relevant to 

(TMTSF)2X compounds. The closing of charge and spin gaps are accurately 
i-^ ■ resolved as topological transitions (jumps in vr of Berry phases). The metal- 

^ , lie or insulating character of each thermodynamic phase is obtained from the 

ground-state expectation value of a displacement operator in reciprocal space. 
PACS numbers: ll.lO.-hx, 3.65.Bz, 7 1.27. -ha 
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Superconductivity studies in the Bechgaard salts with the generic for- 
mula (TMTSF)2X have enriched our knowledge of low-dimensional conduc- 
tors with the new physics discovered.u Although at very low temperatures, 
the transverse hopping should play a role, the extremely large anisotropy 
suggests that a one-dimensional (ID) model should be a good starting point 
for understanding these compounds. Tight-binding calculations suggest that 
only one band crosses the Fermi level.a Therefore, one can construct an ef- 
fective ID one-band model, retaining only the ground state (GS) of each 
unit cell for zero, one, and two holes and can map these states into the cor- 
responding ones of a Hubbard model. This is a formidable task, in contrast 
to cuprate superconductors,^ and has not been done yet. However, it is ev- 
ident that the effective hopping should depend upon the occupation of the 
sites involved, and one expects an important on-site interaction U . These 
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physical considerations lead to the following effective Hamiltonian 

H = Yl (4aCja + h.C.){tAA{l - ■nia){l - Uj^) + tBB rii^Uj^ 

1 1 

+tAB [nia{l - rij^) + nj^{l - Ui^)]} + C/^(nj| - -)(ni| - -) , (1) 

i 

in standard notation, with one hole per unit cell (1/2-filling, n = 1). We 
restrict the present study to the electron-hole symmetric case tAA = tsB = 1- 

Another possible starting point is to retain the states of lowest energy of 
only one organic molecule TMTSF (instead of the unit cell). The resulting 
Hamiltonian is similar to (ll]), but now the operators c\^, c^^ act on low- 
energy states of one molecule (labeled by i), the parameters are different, 
and the appropriate filling is one hole each two sites (1/4-filling, n = 1/2). 

In this work we determine the quantum phase diagram of H, both 
for 1/2- aucl 1/4-filling, using recent developments based on geometrical 
concepts,Ql3'l3 Lanczos diagonalization, and density-matrix renormalization 
group (DMRG) methods. We also extend away from 1/2-filling previous 
studies using continuum limit field theory (CLFT),I3 valid for small values 
of the interactions. Varying the parameters of the model, we obtain several 
thermodynamic phases with noticeable correspondence to the observed ones. 

The basic tool for determining accurate phase boundaries are the charge 
(7c) and spin (7^) Berry phases. 7c (7^) is the phase captured by the GS of 
a system of length L, when twisted boundary conditions cL ^g. = e**'^q^ are 
varied in the cycle < <1> < 27r, keeping $ = cj)^ = <i)j^ ($ = -(^^ = (^^)q 
For systems with inversion symmetry, 7c and 7^ can only be either or vr 
(mod 27r). Then, there are only four possible values of the topological vector 
7 = (70 7s)- In this model, for n = 1, these values correspond to the four 
thermodynamic phases which can be obtained by varying the parameters. 
A jump by vr in 7c (7s) signals the opening of a charge Ac (spin A^) gap 
and a change by e/2 in the polarization (difference between up and down 
polarizations) of the system. For the present H and any filling 7^ = vr if and 
only if As=0, and for n = 1, 7c = if and only if Ac=0. 

For a system of arbitrary length L, in contrast to any other observable, 
7 jumps at the transitions, and the location of this jump in parameter space 
has a very small size dependence. Close to n = 1, the L ^ 00 phase 
boundaries (extrapolated from L = 6, 8, 10, and 12), are displayed in Fig. y 
(a). We estimate the error in the critical value of U to be less than 1%. 

Dotted lines correspond to results obtained using the CLFTJZI Agree- 
ment with the topological transitions is excellent for small t^B — 1 and f7, 
in particular, for 0.6 < Iab ^ 1-2. For large negative U and n = 1, the dom- 
inant correlation functions (CF) at large distances are those for the singlet 
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Fig. 1. (a) Quantum phase diagram for n — > 1. The vector Berry phase 
7 = (7c, 7s) of each thermodynamic phase is indicated. Dotted hues are the 
phase boundaries according to the CLFT. (b) Locahzation parameter |z2,| 
for tjiB = 0.6 as a function of U. Sohd circles are results for L = 12, and the 
solid line is a polynomial extrapolation in 1/L using L = 4, 6, 8, 10 and 12. 



superconductor (SS) and a charge density wave (CDW), as in the negative-C/ 
Hubbard model. For n 7^ 1, SS CF prevail. For < U < 4, n = 1, and low- 
ering tAB, the nature of the different thermodynamic phases is the following: 
a dimerized spin (and pseudospin) Peierls (SP) phase, one with dominating 
spin density wave (SDW) CF, and one with identical triplet superconducting 
(TS) and bond SDW (BSDW) CfE We have established the dominant CF 
of this phase using DMRG, since, by symmetry, the correlation exponent 
Kp = 1 and all CF decay with the same leading power law factor. Slightly 
away from n = 1, where Umklapp processes in the CLFT become irrelevant, 
Kp > 1 and then, the TS CF are larger than all others at large distances. 
Therefore, in materials with weakly coupled chains and n = 1, fluctuations 
in the density of each chain should also favor the TS CF. 

Increasing U, there is a Mott transition TS^SDW which can be studied 
using the quantity zl = {gL\0\gL), where O = exp(i^ Sjo--? "jwl^'^'i \9l) is 
the GS of the system {{gL\gL)=l) with L sites and L particles .BH Note that 
O shifts the wave vector of each particle in 2tt/L. Then, if {gi) describes a 
Fermi or Luttinger liquid with a well-defined Fermi surface in the thermody- 
namic limit, lim/,^00 zl = 0.lj Instead, for an insulator, \zl\ — 1 goes to zero 
in the same way as the Drude weight DcB Furthermore, from a cumulant ex- 
pansion of zl, it can be shown that In |zi| is a measure of the fluctuations of 
the polarization (specifically 2(^)2 In |zl| = (EjaJ^i^)^ - ((Eja J '^jv)^)). 
which are qualitatively different in an insulator and in a metal. The above 
discussion shows that except for some pathological cases, lim^^^oo \zl\ = 
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2. (a) Quantum phase diagram for n = 1/2. Circles and diamonds 
denote the closing of A^ and points where Kp = 1, respectively. Dotted 
lines are the phase boundaries according to the CLFT. (b) Modulus of the 
dominant CF times d? as functions of distance (d) for tAB = 0.6, U = —1. 



in a metal, while lini^^oo |^;l| = 1 in an insulator. In Fig. Mb) we show 
\zl\ as a function of U for a case in which the topological transition is at 
Uc = 2.051. For U < Uc, \zl\ decreases with L, while above Uc {U ~ 3.5 — 4), 
the extrapolated values saturate near 1. 

We studied the evolution of the phase diagram with doping, using nu- 
merical calculations and the extension of the CLFT approachQ away from 
n = 1. For n = 2/3, 1/2, the topological transitions do not completely 
determine the phase diagram, and the calculation of the appropriate gener- 
alization of 7c presents some technical problems. The jump in 7^ still very 
accurately determines the opening of A^. For n = 1/2, the transition to 
dominating superconducting CF at large distances {Kp > 1) is determined 
using numerical data for 8, 12 and 16 sites and finite-size scaling. While at 
n = 1 the computation of Kp has large finite-size effects, for n = 1/2 we 
compared Kp with available Bethe-ansatz results for Iab = IjS with U = —2 
and U = —4 and found the difference to be less than 1%. 

As the system is doped away from 1/2-filling, the SDW advances grad- 
ually over the TS, and a strong increase in the SS region at the expense 
of the former SP takes place. The numerical phase diagram for n = 1/2 
is shown in Fig. ^(0,). Within the CLFT we obtained simple analyti- 
cal expressions Us = 2v2(tAB — 1) for the transition at Kp = 1, and 
Ug = \/2(8/7r + 2){tAB — 1) for the value of U at which A^ closes. As shown 
in Fig. ^(a), these expressions are again in excellent agreement with the 
numerical results, for small U and Iab — 1- However, at intermediate values 
of the interactions, the deviations are much stronger than at 1/2-filling, and 
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the regions in which superconducting CF dominate are larger than those 
predicted by the CLFT. The extension of the SS region for Iab > 1 and 
repulsive U is particularly noticeable. 

The nature of each thermodynamic phase, particularly the TS one, was 
confirmed by direct calculation of CF using DMRG. An example is shown 
in Fig. ^(b). The CF are defined as in Ref. [4]. The slower decay of the TS 
CF for tAB = 0.6, and U = —1 is consistent with the fact that this point lies 
inside the TS phase in Fig. ^(a). 

In summary, we determined the phase diagram of a model, expected 
to describe the low-energy physics of some organic superconductors, using 
recently developed geometrical concepts, and other numerical and analytical 
techniques. The phase diagram is very rich. At 1/2-filling (n = 1), there are 
all possible combinations of open or closed charge and spin gaps, and two 
different Mott transitions take place as a function of the Coulomb on-site 
interaction U. There is a wide region in which triplet superconducting CF 
dominate. This region includes positive values of U if 0.57 < n < 1.43. 

The relevance of the model to (TMTSF)2X compounds cannot be 
demonstrated until the low-energy reduction is performed. However, it is 
remarkable that for n ~ 1 and positive U, the sequence of phases predicted 
by the model with decreasing t^B is the same as that observed with increas- 
ing applied pressure.EJ We expect that n ~ 1 might be more appropriate 
than n ~ 1/2 if the observed dimerization is physically important. 
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